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Abstract 



We report on the results of the ^He('^He,2p)^He experiment at the under- 
ground accelerator facility LUNA (Gran Sasso). For the first time the lowest 
projectile energies utilized for the cross section measurement correspond to 
energies below the center of the solar Gamow peak {Eq=22 keV). The data 
provide no evidence for the existence of a hypothetical resonance in the energy 
range investigated. Although no extrapolation is needed anymore (except for 
energies at the low-energy tail of the Gamow peak), the data must be cor- 
rected for the effects of electron screening, clearly observed the first time for 
the ^He(^He,2p)^He reaction. The effects are however larger than expected 
and not understood, leading presently to the largest uncertainty on the quoted 
5b (-Eq) value for bare nuclides (=5.40 MeV b). 
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I. INTRODUCTION 



Accurate knowledge of thermonuclear reaction rates is important |Q in understanding 
the generation of energy, the luminosity of neutrinos, and the synthesis of elements in stars. 
Due to the Coulomb barrier (height E^) of the entrance channel, the reaction cross section 
(T(i?) drops nearly exponentially with decreasing energy E. Thus it becomes increasingly 
difficult to measure <j{E) and to deduce the astrophysical S{E) factor defined by the equation 



a{E) = ^expi-2nr^), (1) 

with the Sommerfeld parameter given by 2 7r?7 = 31.29 Zi Z2(/i/-E)^/^. The quantities 
and Z2 are the nuclear charges of the interacting particles in the entrance channel, is 
the reduced mass (in units of amu), and E is the center-of-mass energy (in units of keV). 
Although experimental techniques have improved significantly over the years to extend 
(j{E) measurements to lower energies, it has not yet been possible to measure (t{E) within 
the thermal energy region in stars. This region is determined by the Gamow energy window 
Eq ± 6Eq (the Gamow peak) for a given stellar temperature and lies far below the height of 
the Coulomb barrier, approximately at Eq/E^ = 0.01. Instead, the observed (t{E) data at 
higher energies had to be extrapolated to thermal energies. As always in physics, such an 
extrapolation into the unknown can lead to considerable uncertainties. 

The low-energy studies of thermonuclear reactions in a laboratory at the earth's surface 
are hampered predominantly by the effects of cosmic rays in the detectors. Passive shielding 
around the detectors provides a reduction of gammas and neutrons from the environment, 
but it produces at the same time an increase of gammas and neutrons due to the cosmic- 
ray interactions in the shield itself. A An active shielding can only partially reduce the 
problem of cosmic-ray background. An excellent solution is to install an accelerator facility 
in a laboratory deep underground As a pilot project, a 50 kV accelerator facility has 
been installed in the Laboratori Nazionali del Gran Sasso (LNGS), where the fiux 

of cosmic-ray muons is reduced by a factor 10^ 0. The LUNA pilot project (LUNA = 
Laboratory for Underground Nuclear Astrophysics) was designed primarily for a renewed 
study of the ^He(^He,2p)^He reaction {Q = 12.86 MeV) in the energy range of the solar 
Gamow peak {Eq ± 6Eq = 21.9 ± 6.2 keV) for a central star temperature of T = 15.5 x 10^ 
K. The reaction is a member of the hydrogen burning proton-proton (pp) chain 0, which is 
predominantly responsible for the energy generation and neutrino luminosity of the sun. 
So far, the reaction has been studied down to about 25 keV (section II) but there remains 
the possibility of a narrow resonance at lower energies. 

The hypothesis of a low energy resonance was first advanced 0, |10[ as a solution to the 



solar neutrino puzzle, which in those days was regarded as a deficit of B neutrinos. For 



this purpose a resonance with -Er~ 21 keV and F < 6 keV was considered [|11[] as the most 
favourable case. Experimental data available at that time were not inconsistent with the 
existence of a resonance with i?R=15-20 keV and F < 2 keV As more data on solar 



neutrinos became available, it became clear that the deficit of ^Be neutrinos is stronger than 
that of ^B neutrinos. It was shown |T^ that such a pattern of suppression occurs if -Er~ 
21 keV. 
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Such a resonance level in ^Be has been sought |0] without success by various indirect 
routes, and it is also not predicted by most nuclear-structure theories. However, the existence 
of this hypothetical resonance can be positively dismissed only by direct measurements at 
the required low energies (i.e. within the solar Gamow peak). 

For nuclear reactions studied in the laboratory, the target nuclei and the projectiles are 
usually in the form of neutral atoms/molecules and ions, respectively. The electron clouds 
surrounding the interacting nuclides act as a screening potential: the projectile effectively 
sees a reduced Coulomb barrier. This in turn leads to a higher cross section, as{E), than 
would be the case for bare nuclei, cTb(-E), with an exponential enhancement factor []TB|,|T^ 

/lab(^) = a,iE)/a^{E) ~ exp{nr]UjE), (2) 

where Ue is the electron-screening potential energy (e.g. Uc — Zi ■ Z2 ■ e^//2a approximately, 
with an atomic radius). It should be pointed out that for a stellar plasma the value of 
a\,{E) must be known because the screening in the plasma can be quite different from that 
in laboratory studies |]T3|, and <Jh{E) must be explicitly included in each situation. Thus, a 
good understanding of electron-screening effects is needed to arrive at reliable <Jb{E) data at 
low energies. Low-energy studies of several fusion reactions involving light nuclides showed 
H,[TO,IT^ indeed the exponential enhancement of the cross section at low energies. The 



observed enhancement (i.e. the value of Uc) was - in all cases - close to or higher than the 
adiabatic limit derived from atomic-physics models. An exception are the previous ^He+^He 
data (section II), which show apparently a flat S{E) curve down to = 25 keV, although 
the effects of electron screening should have enhanced the data at 25 keV by about a factor 
1.2 for the adiabatic limit Ue = 240 eV. Thus, improved low-energy data are particularly 
desirable for this reaction. 

We report here on such new data obtained by the LUNA Collaboration within the solar 
Gamow peak. Preliminary results, which have been published |TB|, are superseded by the 
present report. 



II. THE 3He(3He,2p)'^He REACTION 

The ^He(^He,2p)^He reaction represents in the exit channel a 3-body-breakup: if the 
breakup is direct, one should observe a continuous energy distribution of the ejectiles de- 
scribed by phase-space considerations; if the breakup follows a sequential process, the ener- 



gies of the ejectiles are described by 2-body kinematics. Experiments have shown P9H21 
that at energies below E = 1 MeV the reaction proceeds predominantly via a direct mech- 
anism and that the angular distributions approach isotropy with decreasing energy. The 
S{E) energy dependence observed by various groups |[T9H23| represents a consistent picture 
(Fig. |1]). The only exception is the earliest experiment |Q where S{E) is lower by a factor 
3 to 5 compared to the other experiments; the discrepancy is most hkely caused by target 
problems ( '^He trapped in an Al foil). 

The absolute S{E) values of refs. pOf and (as well as those from the present work, 
section VI) agree - at overlapping energy regimes - within experimental uncertainties, while 
those of refs. p2[ and are lower by about 25%, suggesting a renormalization of their 



absolute scales. However, in view of the relatively few data points and their relatively large 
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25 to omit 



uncertainties - in comparison to the other data sets - it has been suggested 
these data, without significant loss of information. We verified that no changes in the S{E) 
fit are appreciable (within 1%) by including or excluding the data sets from and p3| . 
A reaction mechanism was suggested at low energies, in which a neutron tunnels 



from one '^He to the other, unimpeded by the Coulomb barrier, up to a radial distance where 
the nuclei overlap appreciably. In this model, a diproton remains and subsequently fissions 
into 2 protons. The calculated energy dependence of the S{E) factor described well the data 
(dotted curve in Fig. p, thus providing confidence in the extrapolation using a polynomial 
function (solid curve in Fig. ffl): 



S{E) = S{0) + S'{0)E + 0.5S"{0)E^ = 5.3 - 3.7 E + 1.9 E^{MeV b) 



(3) 



III. THE LUNA FACILITY 

Technical details of the LUNA setup have been reported Briefiy, the 50 kV accel- 
erator facility (Fig. ^ consisted of a duoplasmatron ion source, an extraction/acceleration 
system, a double- focusing 90° analysing magnet (with adjustable pole faces), a windowless 
gas-target system, and a beam calorimeter. 

The energy spread of the ion source was less than 20 eV, the plasma potential energy 
deviated by less than 10 eV from the voltage applied to the anode, and the emittance of the 
source was 2 cm rad eVi/2^ The ion source provided a stable beam current of about 1 mA 
over periods of up to 4 weeks. 

The high voltage (HV) of the accelerator was provided by a power supply, which has a 
typical ripple of 5 x 10~^, a longterm stability of better than 1 x 10"'^, and a temperature 
coefficient of better than 1.5 x 10^^ /K. The air-conditioned laboratory was kept at a tem- 
perature of T = 21° C and a relative humidity of H = 30%. The HV of the accelerator 
was measured with a resistor chain, contained in an air-tight plexiglas tube, and a digital 
multimeter. The resistor chain was built as a voltage divider, with fifty 20 MQ resistors and 
one 100 kQ resistor (temperature coefficient = 1 x 10~^/K). The multimeter (with a long 
term stability of 5 x 10^^ per year) provided the numerical value of the HV measured across 
the 100 kfl resistor. This HV-measuring device was calibrated at the PTE in Braunschweig 
(Germany) at T = (20 ± 1)°C and H = {35 ± 10)% to a precision of 5 x 10"^ 

The beam entered the target chamber of the differentially pumped gas-target system 
(3 pumping stages) through apertures of high gas flow impedance (C to A in Fig. ||) and 
was stopped in the beam calorimeter. The gas pressure in the target chamber, po, was 
measured with a Baratron capacitance manometer to an accuracy of better than 1%. This 
measurement was absolute and independent of the type of gas used. For ^He gas (99.9% 
enriched in ^He ) of Po = 0.50 mbar pressure, the system reduced the pressure to 1 x 10~^, 
1 X 10~^, and 1 x 10~^ mbar in the three upstream pumping stages; a similar pressure 
reduction was observed for other po values. The gas composition in the target chamber was 
monitored with a mass spectrometer. The ^He gas was compressed by Roots blowers, cleaned 
efficiently using a zeolite adsorption trap (cooled to liquid nitrogen temperature), and fed 
back into the target chamber (gas recirculation). The pressure Po was kept at a constant 
value using a needle valve in combination with an electronic regulation unit. As noted above. 
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the main pressure drop occurred across the entrance aperture A (7 mm diameter, 40 mm 
length, 230 mm distance from aperture B). It was shown that the pressure in the target 
chamber was essentially unmodified by the gas fiow through the entrance aperture A; thus, 
the geometrically extended target zone was characterized by a nearly static pressure. Beam- 



heating effects on the gas density are expected to be less than 0.5% for a maximum 
Po = 0.50 mbar ^He target pressure and the 500 fiA maximum ^He"*" beam current. 

The beam current in the target area was determined to an accuracy of 3% using the beam 
calorimeter (with a constant temperature gradient). The calorimeter was placed (Fig. ^ at 
such a distance d from the center of aperture A [d = 32.2 ± 0.1 cm) that angle straggling 
of the incident beam in the gas resulted in a beam profile smaller than the 200 mm^ active 
area of the calorimeter. 

The LUNA facility was equipped with an interlock system, which allowed the system 
to run without an operator on site. The duty time of the facility in the chosen running 
conditions was about 90%, with a weekly service time of 8 h. As the typical beam current in 
the target area was about 400 /xA, a weekly charge of about 200 Cb could be accumulated 
on the target. 



IV. THE DETECTION SETUP 

The detection setup for the ^He('^He,2p)^He studies had to fulfill the following require- 
ments: 

1. a high absolute efficiency, in view of the expected reaction rates (yields) of about 1 
event/day and less; 

2. a high rejection of natural radioactivity in the detectors, in the target chamber facing 
the detectors, and from the surrounding rocks at LNGS (mainly gammas); tests at 
LNGS have shown that commercial Si detectors exhibited an intrinsic radioactivity 
level, which was about 200 times higher than the above reaction yield; 

3. a high rejection of electronic noise, in view of the needed running times of several 
weeks per energy point; 

4. a clear separation of the reaction products from those of '^He(d,p)'*He {Q = 18.35 
MeV), due to deuterium contamination in the ^He beam (as HD"*" molecules of mass 
3) and in the gas target (found to be smaller). This contaminant reaction has a cross 
section one millionfold higher than that of ^He(^He,2p)^He at Ei^h = 40 keV, mainly 
due to the barrier ratio (d+ ^He)/£'c (^He + ^He) = 0.56, and thus extremely small 
deuterium contaminations (of order 10~^) can lead to sizeable event rates. 

In order to optimize the detection setup and to understand the resulting spectra for quanti- 
tative analyses, a Monte Carlo program was written to simulate the experiment under 



realistic conditions. The Monte Carlo program produces energy and time spectra of the 
ejectiles as well as absolute yields, which could be compared directly with data. Various 
quantitative tests of the Monte Carlo predictions have been carried out successfully p|,|5|,p5 



A Monte Carlo simulation of the proton spectrum from both reactions at a beam energy 
of 50 keV is shown in Fig. ^ a Ni foil of ^He and ^He stopping thickness (20 fim) was placed 
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in front of a 1000 fim thick Si detector (maximum commercially available thickness at the 
time of the experimental project: the desirable thickness was 1400 /im to stop the 14.7 MeV 
protons from ^He(d,p)^He ); a deuterium contamination of 10~^ was assumed. One sees the 
continuous energy distribution of the protons from '^He(^He,2p)^He and the peak(s) plus 
low-energy tail of protons from ^He(d,p)'^He (due to incomplete stopping in the detector). 
Both proton spectra overlap appreciably leading to large uncertainties in the yield analysis 
of the ^He(^He,2p)^He events; thus, the single spectra of Si detectors did not represent a 
viable solution. In addition, singles spectra would not reject the background due to nat- 
ural radioactivity and electronic noise. By the requirement of proton-proton coincidences 
between detectors - placed on opposite sides of the beam axis - a unique signature of the 
^He('^He,2p)^He events could be obtained |^ in the coincidence spectra; also the back- 
ground events (radioactivity and electronic noise) could essentially be eliminated. However, 
the price for this unique signature is a reduction of the absolute efficiency, by at least one 
order of magnitude compared to analysis via single spectra. Therefore, this approach was 
not followed; instead the present approach utilized AE-E detectors, in which single proton 
signals were detected in coincidence between the AE and E detectors of the telescope. In the 
detection setup (Fig. designed according to the indications of Monte Carlo simulations, 
four AE-E telescopes (placed in a rectangular target chamber) were arranged around the 
beam axis: two opposite telescopes each at a distance of 2.7 cm from the beam axis and the 
two other telescopes each at 3.7 cm; the distance from the center of the entrance aperture A 
to the center in front of the telescopes was (8.3 ± 0.1) cm. The different distance from the 
center of the up-down and left-right couples of telescopes was due to mechanical constraints. 
Each telescope consisted of transmission surface barrier silicon detectors with a 0.25 fim 
thick Al layer deposited on both sides of the detectors. The AE and E detectors both had 
an active square area of 2500 mm^; the AE (and E) detector had a thickness of 140 /im (and 
1000 fim) and an energy resolution of 105 keV (and 55 keV) at Ea = 5.5 MeV. A mylar foil 
(1.2 yum thick) and an Al foil (1.5 fim thick) were placed in front of each telescope; they 
stopped the intense elastic scattering yield and shielded the detectors from beam induced 
light. This double shielding was proved, in various background runs, to be effective and free 
from pinhole effects. The detectors were maintained permanently at low temperature (about 
-20°C) using a liquid recirculating cooling system. The ^He ejectiles from ^He(^He,2p)^He 
{Ea, = to 4.3 MeV) and 3He(d,p)^He {E^ = 3.7 MeV) were stopped in the AE detectors, 
while the ejected protons from 3He(3He,2p)^He {Ep = to 10.7 MeV) and 3He(d,p)^He {Ep 
= 14.9 MeV) left signals in both the AE and E detectors of a given telescope (coincidence 
requirement). Fig. ^ shows the identification matrix of one A13-E telescope simulated with 
the Monte Carlo program at -Eiab= 50 keV and assuming a deuterium contamination of 10~^: 
the matrix reveals a clear separation of the events from both reactions. 

Standard NIM electronics were used in connection with the four telescopes. The signals 
were handled and stored using a CAMAC multiparametric system, which allowed for on-line 
as well as for play-back data analyses. Signals from two or more detectors were considered 
coincident within a time window of 1 fis, making negligible the rate of random coincidences. 
A pulser was permanently used in all detectors to check for dead time and electronic stability. 
The acquisition system also stored concurrent information on experimental parameters (such 
as ion beam current /charge, accelerator high voltage, and gas pressure in the target chamber) 
via CAMAC scalers. Controls have also been implemented to stop data acquisition if the 
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beam was lost or to reject an event if an error in the data transmission occured. The system 
handles safely a data flux of about 500 events/s. 

In the analysis of the data, the accepted events from ^He('^He,2p)^He had to fulfill the 
following 3 conditions: 

1. the events had to arise from proton- induced AE-E coincidences in a given telescope; the 
coincidence requirement of each telescope essentially eliminated events due to natural 
radioactivity of the detectors themselves and of surrounding materials; 

2. the events had to lie in a closed region of the AE-E plane (Fig. |^), whose borders were 
first deduced by Monte Carlo simulations and then fixed in order to cut the electronic 
noise; 

3. a proton event had to occur in only one given E-detector; events which triggered more 
than one E-detector were rejected; this condition reduced significantly events due to 
the residual electronic noise, but it also rejected some p-p coincidences. 

The chosen AE-E energy region together with the anticoincidence requirement led to an 
absolute detection efficiency of (7.55 ± 0.15)% as determined by the Monte Carlo program. 
In the AE-E region of events from '^He(^He,2p)^He (Fig. |^), no background events were ever 
observed at LNGS since the installation of the equipment (January 1994) during several "no 
beam" and/or "no target" background measurements (up to two months running) as well as 
during a 10 day run with a ^He beam and a ^He target gas. At the Bochum laboratory (earth 
surface) a background rate of about 10 events/day was found with the same setup. With the 
discussed operating conditions the setup sensitivity in terms of cross section values is better 
than 10~^^ b. Although the selected regions in the spectra of the telescopes allow for a clear 
separation of the events from both reactions (Fig. a few protons from the contaminant 
reaction ^He(d,p)^He can hit the detectors near the edges of their active volumes loosing 
only a fraction of their energy and thus leading to a background rate in the AE-E region of 
the ^He(^He,2p)^He events. The probability of such events was investigated by the Monte 
Carlo program as well as by direct measurements using projectiles with Zj A = 0.5 (selected 
by the 90° analysing magnet). The ratio between the background events in the ^He -^He 
region and those in the clearly separated d -'^He region (= monitor) turned out to be (0.40 ± 
0.04)%. Thus a deuterium contamination equal to 10"'' allows cr(E) measurements as low as 
ii^ = 15 keV, with a signal-to-background ratio equal to 4. During the reported experiments 
the deuterium contamination d/^He ranged between 5 x 10^® and 5 x 10"^. 

V. EFFECTIVE BEAM ENERGY AND ABSOLUTE CROSS SECTION 

For the isotropic emission of the protons in ^He('^He,2p)^He (section II), the number of 
counts dN{z) per unit of time arising from a differential length dz of the extended ^He gas 
target is given by the expression 

dN{z) = NfNy,- cr{E{z)) ■ r]{z) ■ dz, (4) 

where A^t is the ^He target density in units of atoms/cm^ (derived from the measured target 
pressure Po, with a precision better than 1% ), Ab is the number of ^He projectiles per 
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unit of time (derived from the beam calorimeter, with a precision of 3% ), and ri{z) is the 
absolute detection efficiency of all four telescopes including geometry and the acceptance 
criteria discussed in section IV. The efficiency also takes into account that two protons are 
produced per reaction. Introducing the stopping power e (i.e. the energy loss per unit 
length), equation (4) can be rewritten in the form 

dN{E) = NfNi,- a{E) ■ r]{E) ■ t{E)-^ ■ dE. (5) 

The total number of counts for the full target length L = (32.2 ±0.1) cm (from the center 
of the aperture A to the beam calorimeter) is then given by 



N = NfN^ 



j^a{E)-r]{E)-e{E)'^ - dE. (6) 



For the case of a thin target, i.e. introducing an effective beam energy E^s corresponding 
to the mean value of the projectile energy distribution in the detection setup (see below), 
one arrives at 

iV = iVt ■ iVb • c^(^eff) ■ V{E) ■ t{E)-' ■ dE, (7) 

where the values for E^s and the integral were derived from the Monte Carlo program. 

At subcoulomb energies a precise knowledge of the effective beam energy associated with 
observed yields is as important as the yield measurements themselves. For this reason, all 



Monte Carlo predictions have been thoroughly tested [^]. Here they are compared with a 
simple "hand" calculation. Let us consider the case of an incident ^He energy of -Eiab = 50.00 
keV. The absolute energy is known to a precision of ± 2.5 eV, its long term stability is ± 5.0 
eV, and the Gaussian energy spread (FWHM) is 20 eV (section III). When the beam passes 
through the gas target system filled with '^He at the standard pressure of po = 0.30 mbar, the 
beam loses energy in the 3 pumping stages and in the target chamber. For the calculation of 



this energy loss, we used stopping-power values given by the TRIM program [B^]. It should 



be pointed out that experimental energy-loss data at the relevant low energies were 



found to be consistent with those from the TRIM program, within the experimental error 
of 10%. We adopted this error in our analyses. With e = (7.0 ± 0.7)10"^^ eV/atom/cm^ at 
Eiah = 50 keV one finds an energy loss of AEi = (2.2 ±0.2) eV in the 3 pumping stages and 
AE2 = (430 ± 43) eV over the (8.3 ± 0.1) cm distance from the center of aperture A to the 
center of the target chamber (detector location): Ai^tot = (432 ± 43) eV. The quoted error 
arises predominantly from the uncertainty in e . This leads to an effective energy at the center 
of the target chamber of i^iab = (49.568 ± 0.043) keV. The telescopes see an effective target 
length of (7.0 ± 0.1) cm (see Fig. |^). The energy loss over this target length (i.e. the target 
thickness) is AE^ = 365 eV with an estimated error of ± 36 eV. For a constant S{E) factor 
the cross section drops by 11.2 % over this target thickness. Assuming a linear decrease in 
cross section over this region 0, the effective energy evaluation gives E'lab = (49.572 ±0.045) 
keV, or a center-of-mass energy E^g = (24.786 ± 0.023) keV. The error on Ees transforms 
into a ± 1.5% uncertainty on a{Ees). These estimates are consistent with Monte Carlo 
calculations simulating precisely the detection setup (Fig. ^ and the acceptance criteria of 
^He-^He events. The effective energy, as deduced from the mean energy of the distribution 
shown in Fig. |, correspond to -Eiab** = (49.568 ± 0.043) keV. The spread of the energy 
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distribution is determined by various effects: a) the energy straggling p, ^ at the center of 
the target, = 125 eV; b) the thermal Doppler width at room temperature, = 51 
eV; c) the energy spread over the target length seen by the detectors, S^Es = 105 eV. The 
total energy spread 6e is therefore 171 eV, with an estimated error of ±20 eV. This is in 
good agreement with the 189 eV spread predicted by the Monte Carlo program (Fig. 
This procedure was applied for all beam energies and gas pressures. It should be noted that 
the 10% systematic uncertainty in the projectiles energy loss in the target induces only a 
0.1 - 0.2 % error on corresponding to a 1.0 - 3.5 % uncertainty in the S{E) values. 

/^From the geometry of the four square telescopes and the target length L = 32.2 cm 
one expects an absolute efficiency over this length of the order of 10%; the Monte Carlo 
simulations yield (7.55 ± 0.15)%. The geometrical setup efficiency was tested using an a- 
source, which was moved along the beam axis. This test was possible, due to mechanical 
constraints, only recently with a new setup (to be used in future work) where eight 1000 fim 
thick detectors were placed in two consecutive boxes (each containing 4 detectors around 
the beam axis) all at the same distance from the beam axis. The results (Fig. ^ are in good 
agreement with the corresponding Monte Carlo simulations (solid curves in Fig. |^). 

To test the reliability of the corrections applied in the data analyses, we performed 
measurements at Eiat = 50.00 keV using different gas pressures: Po = 0.15, 0.30 and 0.50 
mbar; for each pressure the target thickness, the energy loss, and the detection efficiency 
are different. The resulting S{E) values (Table I) are fully compatible with each other, well 
inside the experimental uncertainties. 

During another test period the detection setup was changed: the detector box was moved 
5.0 cm closer to the beam calorimeter extending the total target length from 32.2 cm to 37.2 
cm. The target pressure was fixed at Po = 0.50 mbar. In this setup a significant decrease (of 
the order of 25%) in reaction yield was predicted by the Monte Carlo simulations compared 
to the standard setup and standard pressure {po = 0.30 mbar), due to the lower effective 
energy. The S{E) results (Table 1) are - within the experimental uncertainties - in good 
agreement with values deduced for the standard setup. 



VI. RESULTS 

Table I and Fig. ^ summarize the ^He(^He,2p)'^He results obtained until December 1996 
with the LUNA setup at the 50 kV underground accelerator facility (LNGS). The lowest 
counting rate was 3 events per day ai E = 20.76 keV. The preliminary data |^ have been 
completely re-analysed and have been integrated with the results of other new measurements. 
The data obtained at higher energies (450 kV accelerator in Bochum) with the LUNA setup 
P are also included for completeness. Previous literature data obtained at E = 24.51 to 
1080 keV |2^, ^ are also shown in Fig. ^. The LUNA data have been obtained at energies 
within the solar Gamow peak, i.e. below the 22 keV center of this peak, and represent the 
first measurement of an important fusion cross section at energies near the center of the 
Gamow peak. No evidence of the hypothetical resonance can be found in the covered energy 
range. The observed energy dependence S{E), for bare nuclides Sh{E) and shielded nuclides 
Ss{E) was parametrized using the expressions P,P^ 

S^{E) = 5b(0) + Sl,{0)E + 0.5 Sl!,{0)E' (8) 
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SsiE) = S^iE)exp{7rr]UjE), 



(9) 



where 5'b(0), 5'b(0), 5*^(0), and f/g are fit parameters. The data set shown in Fig. ^ was 
fitted, in the energy region between 20.7 and 1080 keV, using three methods: 

a) first 5'b(0), 5*^(0), ^'^(O) were obtained considering data for E > 100 keV only; then 
Ue was deduced from all the data and the fixed Sh{E) expression; 

b) all the four parameters were allowed to change and the whole data set was fitted; the 
resulting curves for Sh{E) and Ss{E) are shown in Fig. ^ as dashed and solid curves, 
respectively; 

c) Uc was fixed at 240 eV (i.e. the value calculated in the adiabatic limit) while S'b(O), 
S'b(O) and S{^{0) were allowed to vary. 

The results are summarized in Table II. The three methods give compatible Sb{E) values 
while Ue changes significantly: the methods a) and b) give Ue values higher than the adiabatic 
limit (240 eV), consistent with observations in other fusion reactions. It should be noted 
that fitting the LUNA data only with method b), a screening potential of 490± 30 eV is 
obtained. The difference between observed and predicted Uc values is not understood at 
present. Therefore we adopted S^^Eq) from the theoretical adiabatic screening (method c) 
and derived a further error corresponding to the extreme cases [4 = eV (S'b(-Eo) =5.7 
MeV b) and Ue = 432 eV (^b(^o) =5.1 MeV b): 

Sb{Eo) = 5.40 ± 0.05 ± 0.30 ± 0.30) (Me Vb) 

where the first two errors arise from statistical and systematical uncertainties, respectively, 
and the last error from the lack of understanding of electron screening. 

VII. A HYPOTHETICAL RESONANCE AT LOWER ENERGIES 

The difference between the LUNA data in the 20.7-24.8 keV region and S{E) calculated 
in the adiabatic limit (i.e. fixing f/e=240 eV), can be used to estimate upper limits for the 
strength of a hypothetical resonance located at lower energies. For each measured point we 
calculated an "excess" S-factor, Sexc{E) from: 

5'exc(^) = S^eUE) - S^{E), (10) 

where S'meas(E) are the data given in Table I corrected for a 240 eV screening potential 
and Sad comes from the fit c) in Table II discussed above. The excess cross section, integrated 
in the measured energy interval (20.7 - 24.8 keV), gives the integral yield: 

Fexc = (2.0±0.5)10"^^Me1/ -6. (11) 

We may compare this value with the tail of a low-energy resonance integrated over the 
same energy range. For this purpose, we adopted the following procedure: the resonant 
cross section, (Tres(E), was described according to the usual Breit-Wigner expression: 
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aUE) = rr X^uT .(E)r ^(E)^^-^-^^^^^^, (12) 

where A is the De Broghe wavelength, uj is the statistical factor, E^i is the resonance 
energy and F = Fa + Fb is the resonance total width. Due to the high energies of the 
protons in the exit channel, the partial width of the resonance decay, Fb(-E'), was assumed 
to be energy-independent. For the partial width of the entrance channel, Fa(£'), we used 
the expression p|: 

r.AE) = {—r^Pi{E,R^)ef, (13) 

H tin 

where /i = 1.507 a.m.u. is the reduced mass,i?n = 4 fm is the nuclear radius and the 
dimensionless reduced width 0"^ of the nuclear state has an upper limit of unity (^f < 1, 
Wigner limit). The penetrability Pi{E,Rj^) is given by the equation: 

^'^^'^"^ = F,^(i^,K)lGf(E,K)' (1^) 

and was calculated for s-waves (1=0) using the approximations for the Coulomb wave 
functions, Fi and Gi, given in It turns out that Fa <^ 1 keV at E < 25 keV ( for 9f = 

1) and thus F ~ Fb- 

Equations (12), (13) and (14) lead to the expression: 

11 F 
^-(i^) = '^^F,\E,Rn) + Gf{E,R,,) (E - E^f + (F /2)2' ^^^^ 

where the constant A contains the parameters /i, R^ and Oi- Its value can be determined 
by the comparison with the experimental value l^xc^ 

/■24.8 

{E)dE = Fexc (16) 

J20.7 

With this functional relation 6'f is known for fixed values of the resonance energy E-^ 
and of the total width F. It is then possible to evaluate the reaction rate < av >^cs of 
this hypothetical resonance and compare it with the non-resonant rate < av >ssm used in 



res 



Standard Solar Model (SSM) calculations. Fig. |T0| shows the calculated ratio r =< av > 
/ < av >ssM as a function of E-^ for some values of the total resonance width F, assuming a 
central temperature of the sun of 15 10^ K. The dotted curve in Fig. 10 shows the reaction 



rate ratio in the Wigner limit {6f = 1), which was calculated using equations (15) and (16). 

The < (TV >res valucs obtained with this procedure are upper limits for the reaction rate 
due to the hyphotetical resonance. At energies < 9 keV one concludes that the presence 
of a resonance cannot account for even a partial nuclear solution of the solar neutrino puzzle 



T^,|3^. The same conclusion applies in the energy region between 9 and 20 keV for resonance 
widths > 2 keV. Room seems still available for a very narrow resonance in this interval where 
only direct measurements could rigorously dismiss (or confirm) its existence. 
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VIII. CONCLUSIONS 



For the first time the ^He('^He,2p)^He fusion reaction has been studied in an energy region 
deep inside the Gamow peak, i.e. down to 20.7 keV. This goal has been achieved at the 
new facihty for nuclear astrophysics (LUNA) built in the underground laboratory of Gran 
Sasso (Italy). The results obtained so far show that the S{E) energy dependence is consistent 
with the predictions based on an extrapolation from higher energies. The presence of a low 
energy resonance in the ^He(^He,2p)^He reaction, which could have strong effects on the 
"solar neutrino problem" , does not emerge from the new data. In the near future the LUNA 
collaboration will extend the measurements down to E = 17 keV: the foreseen running time 
here is one year. Definite conclusions with respect to the expected solar neutrino fluxes 
have to await the results of these experiments. For this next phase a new detection set-up, 
designed to reduce the background induced by deuterium contamination in the beam and 
in the gas target, has been developed and is now in a testing phase. 

The electron screening effect, which produces an exponential enhancement of S{E) at 
low energies, can be observed in the new ^He(^He,2p)^He data. This is a second important 
result since previous data did not show a clear evidence of the enhancement (Fig. |1|). The 
screening potential seems to be higher than the value predicted by the adiabatic model of 
the interaction between projectile and atomic electron clouds. It should be noted that the 
adiabatic approximation, among other electron screening models, is that giving the highest 
screening potential. A better quantitative determination of the screening potential will 
be possible when the LUNA collaboration will reach the lowest energy scheduled for the 
'^He(^He,2p)^He experiment. Hopefully ,a complete and clear picture of the screening effect 
will be possible when data collected at very low energies and for other fusion reactions are 
available. This is one of the future goals of the LUNA collaboration. 

The excess of the screening potential could also be explained as a tail of a narrow res- 
onance lying in the not measured low-energy region. Upper limits for the strength of such 
a resonance have been calculated assuming the adiabatic limit for the screening potential. 
Unfortunately, these calculations cannot exclude the existence of a narrow resonance (width 
j 2 keV) at energies between 9 and 20 keV with a strength high enough to give a sizeable 
contribution to the ^He(^He,2p)^He reaction rate in our sun. 
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FIGURES 

FIG. 1. Astrophysical S{E) factor of the ^He(^He,2p)^He reaction as obtained in previous 
work [2^, 2^ (section II). The solid curve is a polynomial fit to the data and the dotted curve a 



theoretical calculation [^6[ normalized to S{0) = 5.1 MeV b. 

FIG. 2. Schematic diagram of the 50 kV LUNA accelerator facility (see also Figs. 2 and 3 of 
§)• 

FIG. 3. Schematic diagram of the rectangular target chamber including the detection setup 
(telescopes) and the beam calorimeter; the given lengths and diameters {(p) are in units of mm (see 
also Figs. H and || in ||5|). 

FIG. 4. Monte Carlo simulation of the protons from ^He(^He,2p)'^He and ^He(d,p)^He detected 
in a 1000 nm thick Si detector (covered with a 20 /xm thick Al foil) at E'lab = 50 keV and in the 
setup shown in Fig. |3|. A deuterium contamination in the '^He beam of d/^He = 10~^ was assumed. 

FIG. 5. Monte Carlo simulation of the AE-E identification matrix for the protons from 
^He(^He,2p)^He and ^He(d,p)^He at E'lab = 50 keV in the setup shown in Fig. ^. A deuterium 
contamination in the ^He beam of d/^He = 10~^ was assumed. 

FIG. 6. AE-E identification matrix of one telescope (obtained in the setup of Fig. ^) at i?iab 
= 50 keV and po = 0.30 mbar. The ^He+^He and d+^He selected regions are shown; note the 
beam-induced electronic noise at the left vertical edge of the matrix. 

FIG. 7. Absolute detection efficiency in a new modified setup (to be used in future work) 
with eight E detectors (arranged in two consecutive boxes) as observed with an a-source {Ea = 
5.5 MeV), which was moved along the beam axis (distance (Ia measured from the center of the 
entrance aperture A). All the detectors have the same distance (2.7 cm) from the beam axis. The 
first section (four detectors) of the setup is nearly equivalent to the configuration in Fig. ^. The 
target length seen by one box (7 cm) is the same as in the setup in Fig. ^. The solid curves show 
the results of Monte Carlo simulations. 

FIG. 8. Monte Carlo simulation of the energy distribution of the '^He projectiles (incident 
energy £^iab = 50.00 keV, ^He gas pressure po = 0.30 mbar) leading to detected/accepted events in 
the setup of the 4 telescopes (Fig. ^). The mean energy of the distribution is -Eiab** = (49.568 ± 
0.043)keV and the energy spread is 5E = 189 eV. 

FIG. 9. The S{E) factor of ^He(^He,2p)^He from the previous work (Fig. |l]) and the present 
work (Table 1). The dashed and solid curves represent Sh{E) and Ss(E), respectively. The solar 
Gamow peak is shown in arbitrary units. The upper right corner inset zooms on the underground 
LUNA data. 
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FIG. 10. The ratio < av >res / < crv >ssM of the resonance and non resonance reaction as a 
function of the resonance energy E'r for some values of the total width F. The dashed curve shows 
the upper limit of this ratio considering the Wigner limit (ef = 1). 
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TABLES 



TABLE I. S(E) factor of ^}ie{^Re,2p)^Re 



Energy ^ 


Charge ^ 


Counts ^ 


Background 


S(E) 




ASsys 


(keV) 


(Cb) 








(MeV b) 




91.70 ^ 


0.0373 


16479 


c 


5.15 


0.11 


0.21 


86.51 ^ 


0.0301 


8931 


e 


5.23 


0.11 


0.22 


81.50 ^ 


0.0544 


8378 


e 


5.33 


0.12 


0.22 


76.29 f 


0.209 


22879 


e 


5.32 


0.11 


0.22 


71.22 f 


0.0945 


5012 


e 


5.43 


0.14 


0.26 


66.06 ^ 


0.0615 


2304 


e 


5.43 


0.15 


0.26 


61.06 f 


0.238 


3562 


e 


5.41 


0.14 


0.26 


55.94 f 


0.257 


2251 


e 


5.50 


0.16 


0.29 


50.64 ^ 


0.825 


2870 


e 


5.63 


0.14 


0.31 


45.82 ^ 


0.784 


1087 


e 


6.14 


0.23 


0.39 


24.80 g 


204.2 


128 


14.5 


5.96 


0.62 


0.23 


24.70 g j 


344.0 


424 


47.5 


6.23 


0.37 


0.24 


24.52 g ' 


49.5 


100 


7.7 


7.10 


0.79 


0.31 


24.30 s 


68.5 


119 


10.9 


6.90 


0.72 


0.37 


24.25 s j 


99.8 


93 


5.0 


6.66 


0.74 


0.26 


23.70 s j 


140.4 


96 


4.2 


6.87 


0.74 


0.27 


23.21 g j 


122.9 


73 


8.2 


7.50 


1.02 


0.30 


23.15 g ™ 


32.7 


28 


2.9 


6.82 


1.47 


0.42 


22.82 g 


139.3 


103 


12.5 


7.21 


0.84 


0.39 


22.78 g j 


307.0 


101 


6.2 


5.97 


0.64 


0.24 


22.33 s 


113.9 


57 


3.7 


7.27 


1.05 


0.40 


22.28 g j 


233.5 


59 


7.3 


5.85 


0.89 


0.24 


21.75 g j 


373.9 


77 


2.4 


7.63 


0.91 


0.31 


21.23 g j 


416.4 


60 


6.5 


7.15 


1.06 


0.29 


20.76 s j 


1044.9 


107 


17.1 


6.80 


0.82 


0.28 



^Effective center-of-mass energy derived from the absolute energy of the ion beam and Monte Carlo 
calculations (including the energy loss of the projectiles in the target gas and the effects of the 
extended gas-target and detector geometries). 
^Deduced from the beam calorimeter (3% accuracy). 
■^Observed events in the ^He-^He region (Fig. ^). 

'^Background events in the ^He-^He region equal to 0.40% of the observed counts in the d-^He 
region of the spectra (Fig. which must be subtracted from the counts in column 3. 
'^Negligible background. 

^Data obtained with the LUNA setup at the 450 kV accelerator in Bochum ||5|, with the standard 
setup and target gas pressures ranging from 0.15 to 0.30 mbar. 

SData obtained with the LUNA setup and accelerator facility at LNGS (energy spread = 20 eV); 
the chosen energy steps are of the order of the energy loss of the beam in the extended gas target. 
'^Statistical error (one standard deviation) including counting statistics and apparative variations 
(pressure (1%), beam power (1.3%), and temperature (1%))- 
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^Systematical error (one standard deviation) including uncertainties in pressure (1%), beam power 
(3%), efficiency (2%), beam energy (200 eV cm. at Bochum, negligible at LNGS), and energy loss 
(10%). 

^Target gas pressure equal to 0.30 mbar, standard detection setup position. 
''Target gas pressure equal to 0.15 mbar, standard detection setup position. 
^Target gas pressure equal to 0.50 mbar, standard detection setup position. 
^Target gas pressure equal to 0.50 mbar, modified detection setup position. 
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TABLE II. S'b(-E) factor of ^He(^He,2p)'^He for bare nuclides and electron screening potential 
energy 11^, from different fitting procedures in the 20.7-1080 keV region. A description of the 
adopted procedures is given in the text. 



5b (0) 


5b(0) 


S^{0) 




x' 


method 


MeV b 


b 


b/MeV 


eV 






5.1 ±0.1 


-2.6 ±0.7 


2.0 ± 1.3 


432 ± 29 


0.93 


a 


5.30 ±0.08 


-3.6 ±0.6 


3.8 ± 1.1 


323 ± 51 


0.85 


b 


5.40 ± 0.05 


-4.1 ±0.5 


4.6 ± 1.0 


240 


0.87 


c 
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Polynomial fit 

S(E)=5.3-3.7E+1 .9E' (MeV b) 
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